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LONG-TERM  GOALS 

The  long-range  goal  of  the  proposed  research  is  to  understand  air-sea  interaction  and  coupled 
atmospheric  and  oceanic  boundary  layer  dynamics  at  low  wind  speeds  where  the  dynamic  processes 
are  driven  and/or  strongly  modulated  by  thermal  forcing.  The  low  wind  regime  extends  from  the 
extreme  situation  where  wind  stress  is  negligible  and  thermal  forcing  dominates  up  to  wind  speeds 
where  wave  breaking  and  Langmuir  circulations  are  also  expected  to  play  a  role  in  the  exchange 
processes.  Therefore,  the  CBLAST-LOW  investigators  seek  to  make  observations  over  a  wide  range 
of  environmental  conditions  with  the  intent  of  improving  our  understanding  of  upper  ocean  and  lower 
atmosphere  dynamics  and  of  the  physical  processes  that  determine  both  the  vertical  and  horizontal 
structure  of  the  marine  boundary  layers. 

OBJECTIVES 

The  goal  of  CBLAST-LOW  is  therefore  to  improve  our  understanding  of  the  processes  that  couple  the 
marine  boundary  layers  under  these  conditions  using  observations,  numerical  simulations  and  models. 
The  ultimate  goal  is  to  incorporate  new  and/or  improved  parameterization  of  these  processes  in 
coupled  models  to  improve  marine  forecasts  of  wind,  waves  and  currents.  Lor  example,  the  boundary 
conditions  between  the  ocean  and  atmospheric  models  are  often  provided  by  parameterizations  based 
on  the  bulk  aerodynamic  method.  This  method  is  also  widely  used  to  estimate  the  fluxes  from  time 
series  measurements  over  the  ocean.  This  briefing  describes  ongoing  efforts  to  improve  bulk  formula 
using  data  collected  during  CBLAST-LOW.  Topics  of  interests  include  a  discussion  of  wind-wave- 
swell  interaction  at  low  winds  and  heat  and  moisture  exchange  in  stratified  conditions. 

APPROACH 

To  achieve  some  of  these  objectives,  the  array  component  deployed  a  3-D  mesoscale  array  to 
simultaneously  observe  the  horizontal  and  vertical  structure  of  the  oceanic  surface  boundary  layer 
south  of  the  tower  as  shown  in  Ligure  1.  This  mooring  component  also  conducted  intensive  ship-based 
surveys  during  the  intensive  operating  period  (lOP).  The  ship-based  surveys  were  coordinated  with 
the  two  aircraft-based  efforts  that  investigated  spatial  variability  of  the  atmospheric  boundary  layer  and 
sea  surface  temperature  field.  The  combined  data  sets  will  be  used  in  conjunction  with  the  modeling 
studies  to  seek  answers  to  unresolved  questions  about  how  the  vertical  as  well  as  the  horizontal 
structure  of  the  coupled  boundary  layers  evolve. 
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CBLAST  2003  Offshore  Array 
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Figure  1.  A  diagram  of  the  CBLAST  region  showing  some  of  the  assets  that  were  deployed  during 

the  main  experiment  lOP  in  the  summer  2003. 


The  tower  component  has  deployed  an  Air-Sea  Interaction  Tower  (ASIT)  spanning  the  water  column 
and  the  lower  22-m  of  the  atmosphere  at  a  water  depth  of  15-m  at  the  Martha's  Vineyard  Coastal 
Observatory  (MVCO)  and  shown  in  Figure  1 .  The  37-m  tower  has  been  instrumented  with  velocity, 
temperature,  conductivity,  pressure,  humidity,  solar  radiation,  turbidity,  precipitation  and  wave 
sensors.  The  tower  is  connected  directly  to  shore  using  a  fiber-optic-conductor  cable,  which  provides 
Gbyte  bandwidth  and  kWatts  of  power  to  the  researchers.  The  velocity  and  temperature  arrays  span 
horizontal  and  vertical  scales  of  O  (1-10)  m  to  resolve  vertical  structure  and  to  permit  separation  and 
quantification  of  processes  associated  with  shear-  and  buoyancy-generated  turbulence,  surface  waves, 
and  Langmuir-like  coherent  structures. 

The  lOP  of  the  main  experiment  was  recently  completed  in  August  of  2003  with  some  components 
continuing  into  the  fall.  The  field  work  during  the  lOP  involved  substantial  collaborations  with  Tim 
Stanton  (NPS)  deploying  complementary  sensors  at  the  ASIT;  Larry  Mahrt  and  Dean  Vickers  (OSU), 
Jielun  Sun  (NCAR),  Djamal  Khelif  (UCI),  and  Haf  Jonsson  (CIRPAS)  obtaining  atmospheric 
measurements  of  turbulent  fluxes,  vertical  profiles  and  horizontal  variability  from  the  LongEZ  aircraft 
in  2001  and  the  CIRPAS  Pelican  aircraft  in  2003;  and  Andy  Jessup  (UW)  and  Chris  Zappa  (LDEO) 
obtaining  IR  remote-sensing  measurements.  In  addition,  we  have  had  substantial  collaborations  with 
regional-scale  modeling  groups  at  Rutgers  University  and  NRE-Monterey,  as  well  as  EES 
investigations  by  Eric  Skyllingstad  at  OSU  and  Peter  Sullivan  at  NCAR.  The  regional-scale  models 
are  providing  a  context  for  interpreting  our  measurements,  and  our  measurements  will  provide  a  means 
of  testing  estimates  of  turbulent  fluxes  and  dissipation  rates  calculated  by  these  models.  The  tower 
measurements  of  horizontal  and  vertical  variability  spanned  a  range  of  scales  similar  to  those  resolved 
by  EES  simulations  and  will  permit  a  quantitative  evaluation  of  EES  model  calculations.  The  proposed 
study  will  produce  a  unique  set  of  simultaneous  measurements  of  turbulent  fluxes  and  dissipation  rates 
on  both  sides  of  the  air-sea  interface,  as  well  as  critical  evaluations  and  improvements  of  turbulence 
parameterizations  used  in  atmospheric  and  oceanic  models.  The  mooring  and  ship  survey 
measurements  spanned  a  range  of  scales  required  to  investigate  processes  on  the  mesoscale  and,  in 


combination  with  the  aircraft  measurements,  will  permit  a  quantitative  evaluation  of  the  eoupled 
mesoseale  model  results. 


WORK  COMPLETED 

Detailed  measurements  of  the  vertieal  structure  of  the  upper  ocean  and  lower  atmosphere  were 
successfully  eondueted  from  the  ASIT  during  the  lOP.  The  atmospherie  arrays  on  ASIT  were 
deployed  in  late  June  and  reeovered  in  early  November,  2003.  During  this  period,  direet 
measurements  of  momentum,  heat  and  mass  fluxes  were  measured  at  3-6  levels  on  the  tower.  These 
measurements  where  eomplemented  by  fixed  sensors  and  a  profiling  package  of  sensors  to  eompute 
mean  profiles  of  veloeity,  temperature,  and  humidity.  Additional  measurements  of  the  radiative  fluxes, 
sea  surfaee  temperature,  preeipitation,  and  the  wave  field  were  eolleeted  to  provide  estimates  of  the  net 
heat  flux  to  the  oeean  and  the  signifieant  wave  height  and  period. 

The  subsurfaee  boom  was  deployed  on  the  ASIT  and  instrumented  during  the  seeond  half  of  the  lOP. 
The  sensors  ineluded  a  horizontal  array  of  ADVs  paired  with  thermisters.  These  measurements  were 
used  to  eompute  subsurfaee  stresses  and  heat  fluxes  during  the  lOP  (Trowbridge  et  ah,  2004;  Gerbi  et 
ah,  2006).  To  obtain  these  fluxes,  a  technique  that  relies  on  differeneing  veloeities  obtained  from 
horizontally  separated  ADVs  is  used  to  remove  the  irrotational  motion  of  the  surfaee  waves.  To  our 
knowledge,  this  is  the  first  eomparison  of  eoineident  direet  eovarianee  Reynolds  stresses  and  heat 
fluxes  measured  on  both  sides  of  the  interfaee. 

The  Nobska  eondueted  4  eruises  during  the  lOP  in  a  wide  variety  of  eonditions.  The  Nobska  was 
outfitted  with  a  direet  eovarianee  flux  system  (DCFS),  IR  radiometers  to  measure  the  SST,  and  a  towed 
thermistor  chain  to  measure  upper  oeean  temperature  strueture  at  very  high  vertical  resolution  during 
transeets  in  the  CBLAST  region.  The  Nobska  also  deployed  a  series  of  drifters  to  doeument  the 
trajeetories  and  the  evolution  of  temperature  strueture  within  a  water  mass.  Some  of  the  towed  and 
drifting  array  results  have  been  processed  and  eombined  with  the  DCFS  results  that  elearly  show  that 
the  surfaee  fluxes  are  rapidly  responding  to  the  spatial  variability  in  the  SST  field. 

Basie  proeessing  and  applieation  of  post-deployment  ealibrations  to  our  data  from  the  2003  lOP  is 
complete  and  further  proeessing  and  quality  eontrol  is  ongoing.  Relevant  portions  of  the  data  have 
been  transferred  to  John  Wilkin  (Rutgers)  for  initialization  and  testing  of  the  high  resolution  Regional 
Oeean  Modeling  System  (ROMS),  to  Shouping  Wang  (NRL)  for  eomparison  with  COAMPS  (Coupled 
Oeean/ Atmosphere  Mesoseale  Predietions  Systems),  to  Larry  Marht  (OSU)  for  eomparison  with 
aireraft  measurements,  and  to  Peter  Sullivan  (NCAR)  for  comparison  with  Large  Eddy  Simulations  of 
wind-swell  interaetions. 

Results  based  on  analyses  of  our  data  have  been  presented  at  the  2004  AGU  Ocean  Scienees  Meeting 
(Farrar  et  al,  2004;  Crofoot  et  ah,  2004;  Edson  et  ah,  2004b;  Hristov  et  ah,  2004b;  Mahrt  et  ah,  2004), 
the  2004  AMS  16*  Symposium  on  Boundary  Eayers  and  Turbulenee  (Edson  et  ah,  2004a;  Hristov  et 
ah,  2004a;  Sullivan  et  ah,  2004;  Trowbridge  et  ah,  2004;  Wang  et  al,  2004a,  2004b),  the  2005  Gordon 
Researeh  Conferenee  on  Coastal  Oeean  Cireulation  (Edson,  2005),  2006  AGU  Oeean  Seienees 
Meeting  (Gerbi  et  al,  2006),  and  2006  AMS  27*  Conferenee  on  Hurrieanes  and  Tropieal  Meteorology 
(Edson  et  ah,  2006;  Sullivan  et  ah,  2006;  Sun  et  al,  2006). 


An  overview  of  the  CBLAST-LOW  program  has  been  aeeepted  to  the  Bulletin  of  the  Ameriean 
Meteorologieal  Soeiety  (Edson  et  al,  2006). 


RECENT  RESULTS 

The  analysis  of  atmospherie  measurements  has  foeused  on  evaluation  of  bulk  aerodynamie  formulae 
and  flux-profile  relationships.  Bulk  aerodynamie  formulae  relate  turbulent  fluxes  of  momentum, 
sensible  heat,  and  latent  heat  to  the  Reynolds-averaged  veloeity,  temperature,  and  humidity,  and  are 
expressed  in  terms  of  dimensionless,  empirieal  eoeffieients  (e.g.  Fairall  et  al.  2003),  whieh  may  depend 
on  quantities  sueh  as  wind  speed  and  wave  age.  Flux-profile  relationships  relate  turbulent  fluxes  of 
momentum,  sensible  heat,  and  latent  heat  to  the  vertieal  derivatives  of  the  Reynolds-averaged  veloeity, 
temperature  and  humidity,  and  are  expressed  in  terms  of  dimensionless,  empirieal  funetions,  whieh 
depend  on  the  ratio  of  the  distanee  from  the  boundary  to  the  Monin-Obukhov  (MO)  length  (e.g. 
Businger  1988).  Bulk  aerodynamie  formulae  and  flux-profile  relationships  are  a  eornerstone  of 
numerieal  weather  predietions,  and  their  uneertainty  is  one  of  the  primary  obstaeles  to  aeeurate  marine 
foreeasts  in  low  to  moderate  wind  eonditions.  The  high-quality  direet-eovarianee  measurements  of 
turbulent  fluxes  obtained  during  CBFAST-low  provide  a  unique  opportunity  for  evaluation  of  these 
relationships. 

Heat  Fluxes:  CBFAST-Fow  measurements  indieate  that  the  standard  TOGA-COARE  3.0  bulk 
aerodynamie  formulation  (Fairall  et  al.  2003)  represents  direet-eovarianee  measurements  of  latent  heat 
flux  aeeurately  when  the  latent  heat  flux  is  positive  (eorresponding  to  an  upward  moisture  flux),  but 
poorly  when  the  latent  heat  flux  is  negative  (eorresponding  to  a  downward  moisture  flux).  Similar 
results  were  reported  by  Edson  et  al.  (2000).  Fieutenant  Crofoot,  a  Navy  student  in  the  MIT/WHOI 
edueational  program,  reeently  eompleted  a  ease  study  of  an  eight-day  period  eharaeterized  by  light 
winds,  a  stably  stratified  atmospherie  boundary  layer,  and  swell-dominated  waves.  The  ease  study 
(Crofoot  2004)  shows  that  failure  of  the  bulk  aerodynamie  estimate  oeeurred  when  adveetion  of  warm 
moist  air  over  eooler  water  resulted  in  a  downward  flux  of  moisture  and  fog  formation  (Figure  2). 

The  uneertainty  in  the  transfer  eoeffieients  for  heat  remains  one  of  the  main  obstaeles  to  aeeurate 
numerieal  foreeasts.  In  partieular,  our  initial  investigations  have  shown  signifieant  differenees 
between  direet  eovarianee  and  bulk  fluxes  in  stable  eonditions,  partieularly  when  the  moisture  flux  is 
direeted  downward  (Fig.  3).  These  periods  of  downward  moisture  flux  are  often  assoeiated  with  foggy 
eonditions.  However,  the  CBFAST  data  indieates  that  the  Dalton  numbers  (i.e.,  the  transfer 
eoeffieient  for  latent  heat)  remain  lower  than  the  COARE  algorithm  parameterization  even  after 
removal  of  downward  fluxes  and  foggy  periods  (Fig.  4).  Therefore,  the  CBLAST  algorithm  proposes 
a  neutral  Dalton  number  that  is  25%  lower  than  the  COARE  algorithm  at  low  winds  speed.  On  the 
other  hand,  the  Stanton  numbers  (i.e.,  the  transfer  eoeffieient  for  sensible  heat)  is  in  reasonable 
agreement  with  COARE  (Fig.  5).  This  result  argues  against  the  eommonly  held  assumption  that  the 
neutral  transfer  eoeffieients  for  heat  and  mass  are  equal. 
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Figure  2.  Time  series  of  the  latent  heat  fluxes  and  visual  evidence  for  the  presence  of fog  during 

periods  of  downward  moisture  flux. 


Latent  Heat  Fluxes 


Figure  3.  Comparison  of  bulk  aerodynamic  versus  direct  covariance  latent  heat  fluxes  measured 
from  ASIT.  The  red  dots  indicate  unstable  conditions  (i.e.,  positive  buoyancy  flux)  while  the  blue 
dots  indicate  stable  conditions  (i.e.,  negative  buoyancy  flux). 
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Neutral  Dalton  Numbers 
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Figure  4.  Individual  and  bin-averaged  estimates  of  the  neutral  Dalton  number.  The  black  lines 
labeled  with  C30  represents  the  COARE  3.0  parameterization  from  Fairall  et  al.  (2003),  while  the 
red  line  represents  the  proposed  CBLAST parameterization. 
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Figure  5.  Individual  and  bin-averaged  estimates  of  the  neutral  Stanton  number.  The  black  lines 
labeled  with  C30  represents  the  COARE  3.0  parameterization  from  Fairall  et  al.  (2003),  while  the 
red  line  represents  the  proposed  CBLAST  parameterization. 


Momentum  Fluxes:  The  transfer  eoeffieient  for  momentum,  i.e.,  the  drag  eoeffieient,  eomputed  from 
the  CBLAST-LOW  data  set  show  good  agreement  in  the  mean  with  the  CO  ARE  3.0  algorithm  (Fairall 
et  ah,  2003),  partieularly  between  4  and  12  m/s  as  shown  in  Fig.  6.  However,  there  is  signifieant 
disagreement  between  the  bin-average  data  and  the  parameterization  at  the  lowest  and  highest  wind 
speeds.  At  low  winds,  we  hypothesize  that  wind-swell  interaetion  is  the  leading  eause.  For  example, 
for  eonditions  with  weak  wind  following  faster  moving  swell,  previous  studies  have  indieated  that  the 
wind  stress  may  be  redueed  relative  to  the  bulk  predietion.  These  eonditions  are  known  as  old  seas  and 
are  eommonly  found  over  the  oeean  whenever  non-loeally  generated  waves  propagate  into  a  low-wind 
region  or  whenever  loeal  seas  slowly  deeay  as  a  storm  moves  out  of  the  region.  The  lowest  panel  of 
Fig.  6  plots  the  bin-averaged  results  for  three  subsets  of  the  data  that  were  measured  in  young 
(developing),  mature  (developed)  and  old  (deeaying)  seas.  If  the  COARE  parameterization  is  eorreet 
for  mature  seas  at  all  wind  speeds,  then  the  bin-averaged  results  indicate  that  the  drag  coefficients  of 
the  younger  seas  are  enhanced  while  those  of  the  older  seas  are  suppressed. 

Although  these  results  as  not  eonelusive,  we  have  reason  to  believe  that  swell  impaets  air-sea  exehange 
at  low  winds  based  on  the  CBFAST  FES  studies  eondueted  by  Peter  Sullivan  (NCAR).  The  FES 
elearly  show  that  fast  moving  swell  in  light  winds  ean  have  a  signifieant  effeet  on  the  wind  field  up  to 
heights  of  O(lOOm).  Under  these  eonditions,  the  wave  driven  winds  produee  a  low-level  jet  and  a  rapid 
deeay  of  the  momentum  flux  with  height.  The  ASIT  data  is  now  being  used  to  verily  these  results  and 
investigate  the  vertieal  strueture  of  the  turbulenee  in  the  surfaee  layer  as  a  funetion  of  wave-age  and 
sea-state.  For  example,  previous  studies  have  shown  that  fully  developed  (mature)  seas  have  a  wave 
age  of  approximately  /  U =1.2  (where  is  the  phase  speed  of  the  dominant  waves),  while 

developing  (young)  seas  have  a  smaller  value  and  deeaying  (old)  seas  a  larger  value.  The  bin- averaged 
profiles  all  depart  from  their  MO  similarity  predietions  as  they  approaeh  the  surfaee,  the  oldest  waves 
showing  a  veloeity  surplus  and  the  youngest  indieating  a  veloeity  defieit  (Figure  7).  These  results  are 
qualitatively  similar  to  the  FES  results  of  Sullivan  et  al.  (2004,  2006). 
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Figure  6.  Individual  (top  panel)  and  bin-averaged  (middle  and  lower  panel)  estimates  of  the  neutral 
drag  coefficient.  The  black  lines  labeled  with  C30  represents  the  COARE  3.0  parameterization  from 
Fairall  et  al.  (2003).  The  bin-average  estimates  in  the  middle  panel  are  also  compared  with  an 
average  parameterization  derived  from  Large  and  Pond  (1981)  and  Smith  (1980)  denoted  by  L&PS, 
and  the  HEXOS parameterization  given  by  Smith  et  al.  (1992)  denoted  by  HEX18.  L&PS  and  C30 
were  developed  using  open  ocean  data  sets,  while  HEX18  was  developed  from  data  taken  in  coastal 
waters  where  the  water  depth  was  18-m.  The  lower  panel  provides  bin-average  results  for  subsets  of 
the  data  that  were  measured  in  young,  mature  and  old  seas  as  characterized  by  the  wave  age 

parameter  Cp/U. 
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Figure  7.  A  comparison  of  the  LES  results  (left  panel)  reported  by  Sullivan  et  al.  (2004)  with 
CBLAST  results  (right panel).  The  LES  results  are  normalized  by  the  value  geostrophic  wind  used 
in  the  simulation  while  the  CBLAST  results  are  normalized  by  wind  speed  at  16  m. 
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IMPACT/APPLICATIONS 

The  2003  lOP  component  of  the  CBLAST  field  program  was  successfully  completed  in  October,  2003. 
Data  quality  and  return  have  been  excellent,  and  a  wide  variety  of  conditions  were  sampled,  including 
low-to-moderate  wind  conditions  and  the  passage  of  strong  atmospheric  and  oceanic  fronts  through  the 
study  region.  The  ASIT  and  the  fifteen  moorings  that  were  deployed  provide  a  complete  time  series  of 
the  passage  of  oceanic  fronts  and  other  processes  with  a  spatial  resolution  on  the  order  of  4  km,  and  the 
ship  based  measurements  complement  this  data  by  providing  a  spatial  resolution  of  about  8  m.  In 
conjunction  with  aircraft-based  measurements  and  satellite  data,  the  in  situ  measurements  collected 
during  the  2003  lOP  constitute  an  unprecedented  record  of  the  evolution  of  the  coupled  air-sea 
boundary  layers.  These  measurements  will  facilitate  a  more  complete  understanding  of  the  relative 


roles  of  local  air-sea  interaction  and  other  processes  (e.g.  ocean  fronts  and  advection)  in  influencing 
the  evolution  of  the  coupled  air-sea  boundary  layer  in  low-to-moderate  winds.  Through  ongoing 
collaboration  with  numerical  modeling  groups,  we  anticipate  that  this  data  and  improved 
understanding  of  air-sea  interaction  will  contribute  directly  to  improving  the  skill  of  marine  forecasts. 

TRANSITIONS 

In  addition  to  several  ongoing  ONR  projects,  the  ASIT  is  being  used  by  investigators  funded  by  the 
NSF  and  NASA  to  conduct  their  research.  The  ASIT  has  become  a  component  of  the  MVCO. 

RELATED  PROJECTS 

James  Edson,  in  collaboration  with  Peter  Sullivan  (NCAR)  and  John  Wyngaard  (PSU),  has  used  the 
ASIT  in  an  NSF  and  ONR  jointly  sponsored  program  entitled  Ocean  Horizontal  Array  Turbulence 
Study  (OHATS):  An  Investigation  of  Subfilter-Scale  Fluxes  in  the  Marine  Surface  Layer.  Detailed 
information  about  this  project  is  provided  in  the  ONR  annual  report  submitted  by  Sullivan. 
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